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Abstract 

The Bochum experimental enhancement of the d+d fusion rate in a deuterated 
metal matrix at low incident energies is explained by the quantum broadening of 
the momentum-energy dispersion relation and consequent modification of the high- 
momentum tail of the distribution function from an exponential to a power-law. 



1 Introduction 



Anomalous enhancement of the sub-barrier nuclear fusion reaction d(d,p)t in a 
deuterated metallic matrix has been experimentally observed at energies of the 
incident beam lower than few keV [1,2]. Electron screening is not sufficient to 
explain this enhancement, and other quantitative explanations are missing [3]. 

In fact, the calculated (adiabatic limit) electron screening potential energy 
Uad is 28 eV, but experiments show an exponential enhancement of the cross 
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section at very low energies that would correspond to a screening energy Uex — 
309 ± 12 eV [1] . Similar behaviors are found in gaseous targets. 

These discrepancies between the calculated Uad and the experimentally in- 
ferred have not yet been understood [3] ; these puzzling results could have 
important consequences also for the study of nuclear fusion in astrophysical 
environments [4,5]. In this paper we discuss a possible explanation of this 
enhancement which is based on the quantum- uncertainty dispersive effect be- 
tween energy and momentum that was proposed by Galitskij and Yakimets 
[6] and recently discussed and applied by Starostin et al. [7,8,9]. 



2 Anomalously large electron screening 

Experimental data for the d(d,p)t reaction in a deuterated Tantalium target [1] 
are reported for beam energies in the range 4-20 keV. The target is cooled with 
nitrogen at a temperature of 10 °C, which corresponds to a thermal energy 
kT = 2.44 X 10-2 eV. Data from Ref. [1] are plotted in Fig. 1. 
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Fig. 1. Measured values for (a) a ■ Vrei and (b) the astrophysical factor 
S = cr{Ecrn) • X cxp ( Eq / Ecni ) • Data are compared to the "bare curve" of 
Eq. (1) and the "screened curve" of Eq. (2) for two values of the electron potential: 
the theoretical adiabatic upper limit Ue = 28 eV and a fit to the experimental data 
Ue = 309 eV [1,2]. 
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In panel (a) of Fig. 1 the experimental values of a{Ecm)vreiiEcm) are plotted 
versus the center of mass energy Ecm- Since the incoming particles are not 
rclativistic and the target particles are practically at rest, these data arc plot- 
ted using Vrei{Ecm) = ^J^Ecm/ 1^' and Ecm = Ebeamf^, where 11 = mD/2 is the 
reduced mass). 

The unscreened cross section can be written as: 

^cm \ V ^cm ) 

with Eg = 2fi Z^Z|e^//i^ . The astrophysical factor S{Ecm) should vary slowly 
in this energy range and is linearly approximated as S{Ecm) = 'S'o + SiEcm , 
where So and Si are extrapolated from energies Ecm > 20 keV (see panel (b) 
in Fig. 1); the values reported in Ref. [1] are Sq — 43 keV b and 5"! = 0.53 b. 

The interacting nuclei "feel" a effective potential barrier lower by an amount 
equal the electron screening potential C/g- the resulting screened cross section 
is: 

Cr{Ecm) = CTbiEcm + Ug) ■ 

When Ue <S -Ecm a correction factor /g can be factorized: 

(7(^cm) - fe ■ MEcm) , (2) 

where 

fe{Ecm, Ue) = exp ( TT W ^ • ) . (3) 

The electron screening potential computed in the adiabatic limit Uad consti- 
tutes a theoretical upper limit for Ue'- 

Ue < Uad . 

For d+d reactions, Fig. 1 compares the "bare curve" of Eq. (1) and the 
"screened curve" {Ue = Uad = 28 eV) of Eq. (2) with the experimental data. 
The screened cross section with Ue = Uad = 28 eV underestimates the exper- 
imental data by about an order of magnitude in the few keV energy range; 
the data could be fitted only by using an unphysical electron potential much 



3 



larger than the adiabatic upper hmit: Ue — 309 eV » Uad- 



S, 



exp 



a{Ecm) E^m exp tt 




) 



Equation (2) impUes 



exp — 



fe ■ S{Ecm) 



and, once again, only the screened curve with Ue — 309 eV fits the data. 

In summary, there is experimental evidence of anomalously high values of low- 
energy fusion cross sections that would require electron screening potentials Ug 
one order of magnitude larger than their adiabatic limit, if explained in terms 

of screening. In deuterated metal targets, the effect depends strongly on the 
metal [2]). Values of Ue > Uad are required also for describing experiments 
with gas target, but violations of the requirement Ug < Uad are less strong 
than in metal targets. 

These anomalous values of f/g are substantially unexplained; the screened po- 
tential approach is probably trying to mimic important processes that have 
been disregarded. One should attempt to find alternative mechanisms that 
could reproduce the enhancement of the cross section at low energy. 



3 Thermal corrections and quantum uncertainty 

The parameters of the experiments (temperature and beam energy) are such 
that thermal corrections can be neglected. In fact if we average avrei over a 
Maxwellian distribution of velocity, the temperature correction factor fx can 
be estimated [3]: 

{crVrel)M = fx ' (crVrel)T=0 , 



For d-|-d reactions {ti'^Eg = 986 keV) even at room temperature {kT ~ 
10~^ keV) and at the lowest energies of the beam {Ei,eam ~ 1 keV) the 
Maxwellian thermal factor ~ 1. 

As showed by Galitskij and Yakimets [6] many-body collisions broaden the 
relation between momentum and energy of the particles. Since momentum 




(4) 
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rather than energy determines the scattering amplitude, the reaction cross 
section must be averaged over a momentum distribution that may differ from 
the energy one. The reaction rate for a mono-energetic beam is: 



{(7Vrel) = j d^Pa Vrel Cr{Ecm)f{Pa) (5) 

oo 

fiPa) = / dE,n{E,) S^iE, - SpJ , (6) 



where the center of mass energy Ecm is function of the beam-particle and 
target-particle momentum, and the momentum distribution function of the 
target particles f{pa) depends on thermal distribution of the target-particle 
energies n{Ea), which we take to be Maxwellian, and on the probability that 
a target particle with energy Eg, has momentum — ^flm^ie^^ (dispersion 
relation), ^^{Ea - ep„). 

Quantum effects are responsible for the dispersion relation between energy 
and momentum S^{Ea — e^^) not being a (5-function, but a broader distribu- 
tion. According to Galitskij and Yakimets [6] the relation between energy and 
momentum is a Lorentzian 



with 7 = hucoii — hn Gcou Vcoii, where Ucou and acoii are the coUisional frequency 
and cross section, while n and Vcoii are the colliding particles density and 
velocity. 

Assuming that the collision cross section could be approximated by the bare 
Coulomb cross section acou — tt e^/ e^^ , the resulting 7 is 



TThNape* hEa 



Pa 



where the relations Vcou = y'^Ea/rria and n — Nap/A have been used with p 
the total density and A the (average) atomic number. 

In the asymptotic regime e^^ ^ /cT, relevant for the particles that undergo 
fusion, the Maxwellian contribution (27rma /cT)"^/^ x exp (— ep^/(/cT')) is neg- 
ligible and the distribution of momenta becomes 



r, . 1 hNApeHT 27TkT 



(27rmaA;r)3/2 Aei 



Pa . 
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Since only particles in the high-energy tail of the distribution e^^ ^ kT con- 
tribute to the fusion rate, the quantum effect contribution 

Q. Pa 



is the only important contribution to the rate. 

More in general, the fact that the relation between energy and momentum is 
not a 5-function results in a distribution of momentum f{pa) with a power- law 
asymptotic behavior in spite of the energy distribution n{Ea) being exponen- 
tial. This power-law tail becomes mostly important for reactions whose cross 
sections select high-momentum particles. 

This quantum contribution can be calculated numerically and we are also 
developing useful analytical approximations: we shall report elsewhere the 
derivation of these approximations and their comparison with the exact nu- 
merical integration. 

In the following we give some preliminary results using a parameterization 
that is motivated by the form of Eq. (10) and that can be used to qualita- 
tively estimate the importance of such quantum broadening of the momentum 
distribution. 
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Fig. 2. The rate fq = {crvj-ei) / io'Vrei)T=Oi where {avrei) is given in Eq. (11) for three 
values of Emax = -E'feeam/10, Eheam/^, and Ehgam/^- For comparison the screening 
correction factor /e of Eq. (3) is also plotted for three values of the electron potential 
Ue = 200, 250, 300 eV. 
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If we use Emax < Ebeam as ail upper bound for ep^ in the integral in Eq. (10) 
(the precise value of the low bound is inconsequential) and work in the relevant 
approximation that the target particles have energies lower than the beam 
particles, we obtain the following partial parameterization for the dominant 
contribution 



In Fig. 2 we have plotted the ratio J'q = {cVj-ei) /{o'Vrei)T=o with (aVrei) given 
by Eq. (11) for three different values of E^ax- The function fq shows an evident 
enhancement at low energy starting from the region of 1-2 keV in qualitative 
agreement with experiments. The threshold below which this enhancement 
becomes important depends on E^ax- We are completing a more detailed 
calculation that does not require the introduction of the parameter Emax and 
that, therefore, can better test the relevance of this quantum effect for the 
experimental results. 



4 Conclusion 

The theory of Galitskij and Yakimets predicts that quantum indeterminacy 
broadens the relation between energy and momentum. 

We have performed a preliminary calculation to estimate the effect of this 
broadening on the momentum distribution of deuteron in metals and, there- 
fore, on the cross section of the reaction d(d,p)t. This calculation shows that 
this quantum effect should give an important enhancement of the cross section 
at low energies similar to the one observed in the Bochum experiments. 

A more quantitative comparison between theory and experiments requires 
the use of more sophisticated analytical or numerical analyses [10] and the 
inclusion of the effects of screening both on the fusion cross section and on 
the ion-collision cross section. This further work is being completed and will 
be published in the near future. 
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